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ABSTRACT 

SWIFT J1756. 9-2508 is one of the few accreting millisecond pulsars (AMPs) discovered to date. We 
report here the results of our analysis of its aperiodic X-ray variability, as measured with the Rossi 
X-ray Timing Explorer during the 2007 outburst of the source. We detect strong (~35%) flat-topped 
broadband noise throughout the outburst with low characteristic frequencies (~0.1 Hz). This makes 
SWIFT J1756. 9-2508 similar to the rest of AMPs and to other low luminosity accreting neutron stars 
when they are in their hard states, and enables us to classify this AMP as an atoll source in the 
extreme island state. We also find a hard tail in its energy spectrum extending up to 100 keV, fully 
consistent with such source and state classification. 

Subject headings: binaries: close — pulsars: individual (SWIFT J1756. 9-2508) — stars: neutron — 
X-rays: binaries 



1. INTRODUCTION 

Nine years after the discovery of the first accret- 
ing millisecond pulsar (AMP), the number and va- 
riety of this type of neutron star low-mass X-ray 
binary (NS-LMXB) continues to grow. The idea 
that a neutron star could be spun up by accre- 
tion t orques in the bosom of a LMXB (lAlpar et all 
119821: iBhattacharva fc van den Heuvell ll99lF was con- 
firmed in 1998, when the first millisecond pulses 
were seen from a LMXB (| Wijnands fc van der Klisl 
1998). Besides coherent pulsations, two types of quasi- 
periodic millisecond X-ray variability have been seen 
in AMPs: twin kilohcrtz quasi-periodic oscillations 
kHz QPOs) in SAX J1808.4-3658 and X TE J1807-294 
Wijnands et al.|[2003t iLinares et al.ll2005l ) and burst os- 
cillat ions also in SAX J1808.4 - 3658 and in XTE J1814 - 
338 (jChakrabartv et all l200l IStrohmaver et alJ 12003( 1 . 
Together, these phenomena place important constraints 
on kHz QPO and type I X-ray b urst models, offerin g a 
new insight onto accretion physics (jvan der Klisll2006l for 
a recent review). As far as low frequency variability is 
concerned, the fastest spinning AMP, IGR J00291+5934, 
showed the strongest broadband noise with the low- 
est char acteristic freq uencies observed to date in a NS- 
LMXB (|Linares et al.ll2007tl . 

Even though AMPs have provided new input to better 
understand accretion onto low magnetic field NSs, impor- 
tant questions about their physics remain open. For in- 
stance, if the magnetic field of the neutron star in AMPs 
disrupts and channels the accretion flow, thereby produc- 
ing the observed pulsations, it is not clear yet why the 
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strength or configuration of such magnetic field should 
diff er drastically from tha t of non-pulsating NS-LMXBs 
(see lCumming et alj|200ll for a possible explanation). In 
this context it is interesting to note that in three systems 
millisecond X-ray pulsations have been seen to appear 
in and disappear from the persistent e mission, produc- 
ing p redominant (HETE J1900. 1-245 5: iGallowav et al. 
2007 1 1 intermittent (SAX J1 748.9-2021:lAltamirano et al. 
200811 or very rare (Aql X-l: ICasella et al.ll2008l ) episodes 
of pulsations. This implies that the AMP within them is 
only active or visible during a relatively small fraction of 
the time, which may provide a link with non-pulsating 
NS-LMXBs. 

On Ju ne 7 th , 2007, a new X-ray transient was dis- 
covered (iKrimm et al.ll2007a( l with the burst alert tele- 
scope (BAT) onboard Swift. Subsequent RXTE ob- 
servations revealed that this was the eighth discovered 
AMP, with a pulse frequency of ~182 Hz and an or- 
bital period of ~54 minu t es dMarkwardt et alJ 120071 : 
iMarkwardt fc Krimml l2007t iKrimm et al.ll2007bD . The 
outburst lasted about two w eeks and a pos s ible in frared 
counterpart was identified (jBurderi et al.1 l2007t l. No 
radio pu lsations were found in quiescence at 8.7 an d 
1.4 GHz ((Possenti et alj|2007t IHessels fc Stappersll2007l) . 

In the present work we analyze the aperiodic variabil- 
ity of the source and compare it with other AMPs and 
wit h atoll sources (a low-luminos ity class of NS-LMXB, 
see lHasinger fc van der Klisl fl989L Using both PCA and 
HEXTE data we also measure its 2-200 keV energy spec- 
trum. We are thereby able to classify SWIFT J1756.9- 
2508 as an atoll source and to characterize its accretion 
state. 



2. OBSERVATIONS AND DATA ANALYSIS 

We used all pointed RXTE observations of 
SWIFT J1756.9-2508 (Table [T]). The observations 
started on June 14 , 2007, and were all performed with 
an ~0.3° pointing offset in order to avoid contamination 
from the nearby (~0.9° away) bright source GX 5-1. 
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TABLE 1 



Log of the observations. 


Set 


ObsID 




Datc a 


Detectors 6 


Count Rate c 512- 


a PDS d 
















(c/s) 




A 


93065- 


-01- 


-01- 


-02 


54265.2 


2.0 


111. Oil. 4 


6 




93065- 


-01- 


■ 02. 


■ 00 


54266.0 


2.0 


102.3±3.8 


16 


B 


93065- 


-01- 


-02- 


-01 


54267.0 


2.0 


95.3±0.9 


4 




92050- 


01- 


■ 01- 


■01 


54267.1 


3.0 


136.7±1.0 


6 




92050- 


• 01- 


■ 01- 


■00 


54267.4 


1.7 


75.3±21.1 


24 




92050- 


-01- 


-01- 


-02 


54268.1 


2.5 


109.3±20.9 


12 


C 


92050- 


-01- 


-01- 


■ 03 


54269.0 


2.0 


71.3±1.3 


6 




92050- 


-01- 


-01- 


■ 04 


54269.2 


2.0 


68.8±4.2 




D 


92050- 


-01- 


-01- 


■ 06 


54270.2 


2.0 


46.0±2.2 


10 




92050- 


-01- 


-01- 


-07 


54270.6 


2.4 


62.4±12.6 


20 


E 


92050- 


-01- 


-01- 


-10 


54271.6 


2.0 


49.4±0.9 


3 




92050- 


01- 


■ 01- 


■ 11 


54271.8 


2.0 


50.0±2.2 


4 




92050- 


01- 


■ 01- 


■ 12 


54271.9 


3.0 


67. Oil. 6 


5 




92050- 


01- 


■ 01- 


■09 


54271.9 


1.0 


20.9±0.3 


5 




92050- 


• 01- 


■ 01- 


■08 


54272.4 


1.3 


28.6±11.1 


19 




92050- 


• 01- 


■ 01- 


■13 


54273.0 


2.0 


47.4±1.2 






92050- 


-01- 


-02- 


-00 


54273.5 


1.0 


21.8±0.6 
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"Observation start date in MJD 
'Average number of active detectors 

c Average and standard deviation of the ~2-35 keV PCA count rate, 
including all active detectors and not corrected for background (which 
we estimate from set E to be ~22.3 c/s/PCU; see Section[2} 

d Number of power density spectra extracted from the observation, 
each of them 512 s long 

After applying standard screening criteria 6 we extracted 
energy spectra from PCA (using only PCU2) and 
HEXTE (using cluster B) data and power spectra from 
PCA data, using all active PCUs. 

After checking that no strong source falls within the 
field of view of HEXTE-cluster B 7 , we extracted source 
and background spectra from standard mode data. The 
background spectrum was obtained by averaging both 
rocking positions. We applied dead-time correction and 
created response matrices taking into account the above 
mentioned pointing offset. 

A first estimate of the PCU2 background was given 
by the PCA background model for faint sources. Given 
that the source lies in the direction of the Galactic bulge 
and from the observed flux in the last part of the out- 
burst (or "tail", set E, see below) it is clear that there 
is additional (non-instrumental) background flux, most 
likely originating in the Ga lactic ridge (note tha t the 
source was not detected by iCampana et al.l 120071 with 
the Swift-XRT on June 21 st , within set E). The pres- 
ence of an iron line at ~6.5 keV whose energy, width 
and flux remained constant as the overall flux decreased 
also suggests that the Galactic ridge is the main source 
of contaminating flux. In order to take this into account 
and under the assumption that both instrumental and 
astronomical backgrounds remain constant during the 
observations, we obtained the average spectrum of the 
tail of the outburst and subtracted it from every PCU2 
spectrum. This final empirical background estimation is 
independent of PCA background modelling. The iron 
line was not present in the resulting spectra. We created 
response matrices for the ~0.3° off-axis pointing. 

A detailed modelling of the energy spectra is beyond 
the scope of this paper, and probably difficult due to 

6 See http://heasarc.nasa.gov/docs/xte/abc/screening.html A 
conservative filter on the time after BAA passage (>15 min) and 
the ELECTRON2 rate (<0.1 c/s) was applied before extracting 
the energy spectra. 

According to HEXTErock: 



the faintness of the source, so we opted for the fol- 
lowing phenomenological approach. We fitted the joint 
PCA+HEXTE spectra with an absor bed blackbody plus 
power law model using XSPEC v.ll (jArnaudll!99(l . al- 
lowing for a constant factor (weighted mean 0.68±0.03) 
to account for the difference in effective areas between 
PCA and HEXTE. The equivalent hydrogen column den- 
sit y was fixed t o 5.4xI 22 cm -2 , the value measured 
bv lKrimm et a l. (2007b) from the Swift X-ray telescope 
(XRT) spectrum. Figure Q] displays the evolution of 
the 2-200 keV unabsorbed flux and the BAT 15-50 keV 
lightcurve, showing the rise, decay, duration and peak 
flux of the outburst as well as the data sets used in the 
timing analysis (see below). 

In order to study the low-frequency variability, we per- 
formed fast Fourier transforms (FFTs) on 512 s event 
mode data segments keeping the original 125 /is time 
resolution. The power spectra thereby cover the ~0.02- 
4000 Hz frequency range. We used PCA channels to 
79 (^2.0-35.0 keV) excluding the highest (background- 
dominated) channels to increase the signal-to-noise ratio, 
and averaged all power spectra within each observation. 
No dead time nor background correction was made be- 
fore the FFT s and the Po i sson n oise level w as subtracted 
accor ding to iKlein-WoTtl (|2004f ) (using the I Zhang et~aLl 
119951 formula and shifting the resulting level to match 
the 2000-4000 Hz range, where no intrinsic power is ex- 
pected). When looking for changes in frequency and 
power of the broadband noise (see Secj3|) we averaged 
observations successive in time and with power spectra 
consistent within errors in order to improve the statis- 
tics. The outburst was divided into five contiguous data 
sets (labeled A-E; see Table Q] and Figured]). In set C 
low statistics only allows us to constrain the break fre- 
quency and the flat-top power level. In data set D the 
source was too faint to measure accurately its energy or 
power spectrum (about 2.5 c/s/PCU), and as discussed 
above data set E is backgro und dominated. P ower spec- 
tra were "rms-normalized" (Ivan der Krisf l995) using the 
count rate observed in set E as background count rate. 

We fitted the broadband power spectr a with a sum 
of fou r Lorentzians in v max representation ( jBelloni et al.l 
[2002] ) 8 . three of which were zero-centered. Following pre- 
vious notation dBelloni et all 120021: | van Straaten et aLl 
120031 iLinares et all 120071 ,. we call the different 
Lorentzians Li and refer to their characteristic frequen- 
cies as Vi 7 where the label "i" identifies the variability 
feature (see Section [3]). 

3. RESULTS AND DISCUSSION 

Figure [2] shows the broadband power spectrum of 
SWIFT J1756.9-2508. Its overall shape is typi- 
cal of island and ext reme island state s (IS/EIS) of 
atoll sources (see e.g. iBarret et al.l 2000: B elloni et al.l 
2002), as well as similar to low (luminosity) and hard 
(spectrum) states of black hole X-ray binaries (e.g. 
XTE J1118+48 0: iRevnivtsev et~all [2000. XTE J1550- 
564; iCui et alJ I1999T I. Furthermore, the integrated 



Vmax = y v o + A 2 gives the characteristic frequency of the 
variability feature, where vq is the Lorentzian's centroid frequency 
and A its HWHM. The quality factor Q = vq/2A is a measure of 
the coherence of the feature. Its strength is given by the integral 
power (0— oo) whose square root, in the normalization we use, is 



http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/HEXTErock/HEXTErock.plthe fractional rms amplitude of the variability. 
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(0.01-100 Hz) fractional rms variability was very high 
(~35%), also characteristic of such low- hard (Comp- 
tonized) states. This allows us to classify this AMP as an 
atoll source, so mething commo n to all the other AMPs 
known to date (Wiinands 2005). The low characteristic 
frequencies of the variability (see below) indicate that 
SWIFT J 1756.9-2508 was in the EIS during its 2007 out- 
burst (see Ivan der Klidl20M for an overview of states). 

TABLE 2 

Parameters of the 4-Lorentzian fits to sets A and B. 



Par 


amctcr 


Lh 




L l>n„, 


L u 


X 2 /d.o.f. 


Set 


A 














a (Hz) 


0.12±0.01 


0.45±0.03 


3.1±0.9 


53±30 


167/150 




(%) 


17.7±0.8 


13.1±1.6 


17.9±1.3 


21±3 




Q 




(fixed) 


0.9±0.3 


(fixed) 


(fixed) 




Set 


B 














IX (Hz) 


0.09±0.01 


0.39±0.02 


3.1±0.6 


141±41 


142/150 




(%) 


18.3±0.8 


15.1±1.5 


21.0±0.9 


34±3 




Q 




(fixed) 


0.7±0.2 


(fixed) 


(fixed) 





See Section [2] for details of the fit function and notation. All errors 
are la. 

The most prominent variability feature is a strong, 
band- limited and flat-topped noise component (Lb, 
for "break"), with a break or cutoff at a frequency 
^ o ~0.1 Hz. This break frequency lies in the low frequency 
end of the range usually observed in NS-LMXBs. We also 
detect a narrower (Q~0.8) component at higher frequen- 
cies (~0.4 Hz). The correlations between characteristic 
frequencies of power spe ctral features are a useful tool 
to identify such features (IWiinands fc van der KlisllT999t 
iPsaltis et ai] 119991 Ivan Straaten et all 120051) . The nar- 
ro w component at ^0.4 Hz follows the correlation found 
by IWiinands fc van der KlTi (|1999| ) between Vb and the 
frequency of the "hump" (Lh) present above Lb in both 
black hole and neutron star systems (see Figure [3]), and 
we therefore identify it as Lh- The next component in 
order of increasing frequency has a characteristic fre- 
quency of 3.1 Hz, close to (but a bit lower than) what 
we would expect from extrapolating t he relation between 
Llow and L b found in atoll sources (|van Straaten et ahl 
120021 120031 ). This leads us to tentatively identify it as 
Ltow The highest frequency component that we find 
in the power spectrum of SWIFT J1756. 9-2508 is most 
likely L Ul also present also present in the same frequency 
range in IS/EIS of atoll sources. A less likely candidate 
is the so-called "hecto-Hertz" Lorentzian, as that compo- 
nent seems to be present only when the variability is at 
higher frequen cies, with Vb more than o ne order of mag- 
nitude higher (|van Straaten et afll2002h . All parameters 
from the fits to sets A and B are shown in Table O 

In disk accreting X-ray binaries as the flat-topped noise 
shifts down in frequency its flat top becomes higher. In 
other words, the rms level at the break of the power 
spectrum is anticorrelate d with the frequency of t hat 
break. This was found by iBelloni fc Hasinger] (|1990t ) in 
Cyg X-l and extend ed to other BH and NS systems by 
IBelloni et all (|200^ . SWIFT J1756.9-2508 also follows 
this anticorrelation, as can be clearly seen from Figure [3J 

The broadband 2-200 keV spectra from all observa- 
tions within sets A, B and C were reasonably well fitted 
(reduced \ 2 < 1-5 for 61 degrees of freedom) with the 
adopted (Sec. [2]) model (see Figure [4]). The resulting 



parameters are shown in Table [3l A hard tail extending 
up to ^100 keV was present in the energy spectra during 
most of t he outburst, similar to island states of other atoll 
sources (|di Salvo et al.ll200l . As can be seen from Fig- 
ure [H the resulting power law index increases when the 
outburst decays. Several models have been proposed for 
the emission process at work in the transition from out- 
burst to quiescence, including pulsar sho ck emission and 
therm al emission from the neutron star (jCampana et al.l 
119981 for a review). Our data are not sensitive enough to 
yield precise tests of such models, and in general sensitiv- 
ity limits and time constrains make this task a difficult 
one. 



TABLE 3 
Parameters of the spectral fits. 



ObsID 




Ph. I. 


P.L. 


kT 


x 2 /d.o.f. 


2-200 keV una 


bs. 










(keV) 




Flux (erj/s/e» 


> 2 ) 


Set A 


93065-01-01- 


02 


1.85±0.02 


(13.0±0.6)X10 -2 


0.5±0.2 


43/61 


(13.0±0.4)xl0 _ 


10 


93065-01-02- 


■00 


1.84±0.02 


(11.0±0.6)X10 -2 


0.7±0.1 


53/61 


(11.7±0.2)xl0~ 


10 


Sct B 


93065-01-02- 


01 


1.92±0.04 


(11.0±0.9)X10 -2 


0.5±0.1 


59/61 


(9.7±0.5)X10~ 


10 


92050-01-01- 


■01 


1.94±0.04 


(11. Oil. 1) X 10~ 2 


0.6±0.4 


61/61 


(8.7±0.4)xl0~ 


10 


92050-01-01- 


■00 


1.88±0.02 


(9.6±0.5) X 10~ 2 


O.BiO.l 


70/61 


(8.0±0.2) X 10~ 


10 


92050-01-01- 


■ 02 


1.90±0.03 


(9.1±0.6) X 10 -2 


0.7±0.1 


83/61 


(6.8±0.3)xl0" 


10 


Set C 


92050-01-01- 


■ 03 


2.03±0.05 


(6.4±0.6)X10 — 2 


0.5±0.1 


73/61 


(2.8±0.3)X10" 


10 


92050-01-01- 


■ 04 


1.98±0.06 


(4.7±0.7) X 10~ 2 


0.4±0.2 


41/61 


(2.8±0.3) X 10~ 


10 


NOTE.- 


- In all fits the column density was 


fixed to 5.4xl0 22 cm 


-2 



(Krimm et al. 2007b). The power law (P.L.) is characterized by the 
photon index (Ph. I.) and the normalization, which corresponds to the 
photon flux density at 1 keV. In five cases the B.B. normalization was 
not significant but in all 8 observations it gave a 2 — 4 a improvement 
in the fit, according to the f-test for additional terms. All errors are 
la. 

We do not detect a significant change in the frac- 
tional rms of Lt with energy in the studied energy 
range, suggesting that the spectrum of this component 
is flat. Previous work has shown that in LMXBs the 
X-ray luminosity does not uniquely de termine the char- 
acteristic frequencies of the variability dFord et alJ[2000b 
Ivan Straaten et all 120051 : Ivan der Klisl 12006] ) so that an- 
other parameter must act to decouple luminosity and 
variability frequencies. Figure [5] shows that as the out- 
burst decays and the inferred mass accretion rate drops 
the fractional rms variability seems to increase, whereas 
the frequency of the break decreases. The same be- 
havior has been ob served e.g. in IGR J00291+5934 
(|Linares et al .1120071) . A successful model for broadband 
noise in LMXBs should explain both this short-term cor- 
relation between luminosity and break (or "cutoff") fre- 
quency and the long-term luminosity-frequency decou- 
pling discussed above. 
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Fig. 1.— Light curve of the outburst of SWIFT J1756.9-2508, showing 2-200 keV and 15-50 keV fluxes measured with RXTE- 
PCA&HEXTE and Sm/f-BAT , respectively. Filled circles show fluxes after removing the Galactic ridge contribution (see Section [2] for 
detail s), while for open squ ares only the instrumental background was subtracted. The upward arrow indicates the discover y date (June 7 th , 
2007; Krimm et al. 2007a) and the downward arrow the date of the first non-detection with Sun/t-XRT (June 21 s ', 2007; Campana et al. 
120071) . Data sets A-E used in this work are shown with their corresponding labels along the bottom axis. 
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Fig. 3. — Uppe r panel: Correlation between the characteristic frequencies of and L^, after Wiinands & van der Klisl (1999) ( see 
IBellon i ct al. 200j for a discussion and extension of this corr elation). Lower panel: Co rrelation between the frequency of the break in the 
power spectrum and the power level at the break frequency (Bclloni & Hasingcr 1990). SWIFT J1756. 9-2508 follows both of them (see 
Sec. |3]for details). 
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Fig. 5. — Time evolution of spectral index (top), integrated (0.01-100 Hz) fractional rms (middle) and frequency of the break in the 
power spectrum (bottom). The BAT light curve is shown in grey and the data sets used are indicated along the bottom axis. 



